Levy stable two-sided distributions: exact and explicit densities for asymmetric case 
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We study the one-dimensional Levy stable density distributions g{a,j3;x) for — oo < x < oo, for 
rational values of index a and the asymmetry parameter /3: a = l/k and (3 = [l — 2r)/k, where I, k 
and r are positive integers such that < l/k < 1 for < r < Z and 1 < l/k < 2 for I — k < r < k. We 
treat both symmetric (/? = 0) and asymmetric (/3 7^ 0) cases. We furnish exact and explicit forms of 
g{a, P; x) in terms of known functions for any admissible values of a and /3 specified by a triple of 
integers k, I and r. We reproduce all the previously known exact results and we study analytically 
and graphically many new examples. We point out instances of experimental and statistical data 
that could be described by our solutions. 
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The probability distributions characterizing anoma- 
lous diffusive behaviour fiave been a subject of intense 
activity on experimental and theoretical side. Among 
various forms proposed, that of the heavy-tailed Levy 
stable laws is of widespread use, mostly due to their uni- 
versal presence in such diversified fields as econophysics 
P, Q, physics of amorphous materials Q, geology [1], 
biophysics Q, statistics of phone networks internet 
traffic and dynamics of human contactsl. Entire 
monographs are devoted to the thorough study of this 
huge field Comprehensive reviews are available re- 
porting on the state of the art Il(>jl3j . Further references 
may be traced back from [2l| . 

The goal of the paper is to give exact and explicit 
expression for the two-sided Levy stable distributions 
g{a,l3;x), < a < 2, for symmetric (/3 = 0) and asym- 
metric (/? ^ 0) cases. Since g{l,l3;x) requires special 
treatment [9| we omit the case a = 1. 

The probability density function (PDF), g{a,f3;x), is 
called stable if the product of characteristic functions 
(CF) of two such laws is a CF of another law of the same 
type. The general PDF of this type g{a, /3; x), where ei- 
ther — 00 < a; < cxD, or a; is confined to one of semiaxes, 
see below, has the CF defined as the Fourier transform 
in the form [i,[lM3 



g{a,l3;p) = exp cxp (i/37rs5n(p)/2)] 



(1) 



where \p\ and sgn{p) are the absolute value and the 
sign of p, respectively, and g{a,(3;p) satisfy the rela- 
tion g{a,l3] —p) = g{a, —j3]p). According to the values 
of parameters a and /3 we can distinguish the following 
variants: (i) for < a < 1 and < a: we have for 
(3 = —a one-sided PDFs defined for < a; < oo, whereas 
for 13 — a they are defined only for — cxd < a; < 0, oth- 
erwise they are two-sided, i. e. —00 < a; < 00; (ii) for 
1 < a < 2 and < 2 — a the PDFs are always two-sided. 



Only under these restrictions on a and (3 the positivity 
of g{a, 13] x) for all allowed x is guaranteed [1, Hill. 

The finding of exact and explicit form of two-sided 
g{a, /3; a;) turned out to be a true challenge. In the litera- 
ture we can find only a limited number of exact formulae 
for g{a,l3;x). The two-sided asymmetric cases include 
a = 1/3 and /3 = ±2/3 @, a = 3/2 and ^ = ±1/2 
9I jlSl. The symmetric cases {(3 = 0) concern a — 4/3 
Uliai, a = 5/4 and 6/5 [H, a = 1/3 and 1/2 [HHa], 
a = 1/5, 1/4, 2/5, 3/5, 3/4, and 4/5 [l|. The exact 
solutions for one-sided case [3 — ~a for rational a have 
been recently obtained in [21[ . For any other values of a 
and (3 for two-sided situation the only source of informa- 
tion are numerical calculations, often problematic if not 
impossible for small values of a [22|, [23| ■ 

In what follows we shall indicate how the approach 
of pT| can be extended to obtain new exact represen- 
tations for two-sided case with rational values of a and 
/?. One should be reminded at this point that the func- 
tional structure of g(a. 13' x) depends on an essential way 
on the value of a [3, for 1 < a < 2 g{a,/3;x) is 
a unique function for both signs of x. On the contrary, 
for < a < 1, /3 7^ 0, —a, the function g{a,l3;x) is ob- 
tained by matching of two different functions g{a, —j3; x) 
for a; < with g{a^f3;x) for a; > 0, at the point x = 0. 
In any case it is sufficient to use g{a,f3;x) defined as 



1 f°° 

g{a,l3;x) = - K / dp 
Jo 



(2) 



Before going to the most general case let us embark upon 
an intermediate situation where for a given a only certain 
values of /3 intervene. 

The general form of one-sided g{a, —a; a;) for rational 
a = l/k, I < k, where l,k arc integers, was recently 
presented in [21| where it is denoted by ga(x). As an 
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initial approach to the two-sided case we shall generate 
certain two-sided solutions /?*; a;), a* = a~^,l < 

a* <2 and 2, via the duality law (DL) [1 ^ 

applied to one-sided g{a, ^a; x). The DL implies, for 



a < 1, g{a*, 



■2;x) 



-(l+l/a) 



■ 



'i/a-j from 



2 



where for rational a 

to nn 

gia*,a* ~ 2;x) 
fc-i 



l/k we have, according 



l.A{l,l+jl/k) 
A(fcj + 1) 



(3) 



In Eq. (jSj pFq ([b")!^) is the generahzed hypergeomet- 
ric function with the upper and lower parameter lists 
equal to (cp) and {bq) respectively [l^l, and A(n,a) = 



g+l 



, ,^ , ,^ is a special list of n parameters. The 

numerical coefficients bj{l, k) in Eq. ([3]) have the form of 
Eq. (4) in [2l[ . We emphasize that while the distributions 
g{a, —a; x) are defined for x > 0, g(a*, a* — 2; x) are au- 
tomatically valid for — oo < x < oo. To carry out this 
procedure we consider Eq. (6) of [2l| for p = 4 which cor- 
responds in our notation to g(4/5, —4/5; a;), x > 0. The 
DL yields 5(5/4, — 3/4;x) for —oo < x < oo whose exact 
form is given by Eq. ([3]) for I = 4 and fc = 5. In view of 
previous remarks and the values of a* and /3* involved, 
(7(5/4, —3/4; x) is the unique function for both semiaxes. 
In Fig. [T]wc compare both of these distributions. As far 
as we know the function ^(5/4, — 3/4; x) is a new exact 
solution for two-sided case. According to Eq. ^ it can 
be represented, after appropriate simplifications of their 
parameter lists (op) and [bq), as a sum of four hyperge- 
ometric functions of type 3F3 of argument (-4^x^/5^). 
For a given < a < 1 this DL route based on one-sided 
solutions in [2l| will always yield special two-sided den- 
sities in the form g (i, ', x) , —00 < x < 00. They 
should be hitherto considered as known 123 . 



3; x) 



g(4/5, -4/5; X) 




FIG. 1. Comparison of the one-sided distribution 
g(4/5, —4/5; x) and its dual two-sided counterpart 
g(5/4,-3/4;2;), see Eq. ®. 



In fact; we have achieved a more ambitious goal by 



finding exact 5(0;, /3; x) without restrictions on /? imposed 
by the DL. By extending the method of 121| we have 
established a general and universal formula for g{a, /?; x) 
which encompasses both one- and two-sided cases. The 
general form of Eq. ([2]) for appropriate rational a = l/k 
and /3 = a — 2r/fc, where Z, k and r are positive integers 
(see [2^), is given by the exact expression: 



g(a,/3;x) 



m+l^M ^A(A/j + l)" 



'^' c,il,k,r) 

l/M 

l,A{m,l+j m/M) 



i-iy 



-M ' 



, (4) 



where m = min(/, k) , M = max(/, fc), the lower sign is for 
I < k and the upper sign for I > k, with the coefficients 



Cj {I, k, 



Ml/2-J ^l/2+j7n/M 

nLir(^)] [n 



i=0 
M 



J_ 

M 



r ( '~-^~ 



i + l 



n::^[sinH-.^)]- 



(5) 



Here is the sketch of derivation of Eqs. ^ 
and ([5]). They result from the application of 
the Mcllin transform to g{a,f3;x) of Eq. 
A4 [g{a, /3; x); s] = dx x^~^ g{a, I3;x) for complex 

s, which is equal to ^ cos (ir^^ + ir T{.s) F (i^) . 

Then, g{a,P;x) will be perceived as the inverse 
Mellin transform, i. e. it is formally equal to 



M 



The next 



J_ cos f 7r:f + TT s^) r(s)F (i^) ; ; 

steps involve the use of Euler's reflection formula for cos- 
inus, the passage to rationals a = l/k and (i = [l — 2r)/k 
and the use of Gauss-Lcgendre multiplication formula 
for all gamma functions. Putting all the terms together, 
we employ, as an intermediate tool, the storing of 
the inverse Mellin transform as the Mcijer G function 
G™^" {z\^^^) ^^^^ Thc hual use of conversion formula 
8.2'.2.3, p.' 618 of Ref. ^ yields Eqs. Q and ©. 

The actual construction of g{a^ /3; x) from Eqs. Q and 
([S]) for a given triple (Z, k, r) boils down into three distinct 
alternatives: (a) a = l/k < 1, (3 = a — 2r/k = —a which 
gives r = I; it yields a one-sided g (|, —^',x) for x > 0, 
elaborated in [21|; furthermore, for a = Z/fc < 1, but for 
(3 = a — 2r/k = a, which gives r = 0, it yields a one- 



sided glj:, 



defined only for x < 0; (b) for a = Z/fc 



I = \{l-2r)/k\ < 2 -l/k, 
k; both two-sided density 
are defined on 



such that I < l/k < 2 and \j: 
which implies I — k < r < 
functions g(|,-^;x) andg(|,^;x) 
—00 < X < 00 (and arc mutually symmetric with re- 
spect to (0, y), see Fig. [2]); (c) a — l/k < 1, but 
\/3\ = \(l - 2r)/k\ < a, which implies < ?- < /; 
here the density g{a,/3;x) decomposes into two differ- 
ent functions according to the sign of x, and is given 
by 9 (i, ^) 0{-x)+g (I, 21^; x) 0{x), with 0{x) the 
Heaviside function. The matching at x = of these two 
components assures the continuity of 5(0;, /3;x) at x = 0, 
along with all its higher derivatives, see Figs. [21 El and 21 
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Eqs. (HI and (O can be equivalently represented as 
a single infinite series derived by Bergstrom and Feller 
flEi which is a two-sided variant of the Humbert 
expansion [28| . This formula {vide Eqs. (4) and (6) in 
jlSl p is however very slowly convergent in both regions 
of a — > and a < 1. On the contrary, the Eq. ^ 
is easily adaptable to computer algebra systems, with 
built-in ping's providing improved convergence, see [29| 
for ready-to-use Maple® procedure L2S. In Fig. [3] we 
present g(a, /3; x) for small values of a and /? (a = 1/15 = 
—f3, and /3 = 1/45) for which neither Bergstrom- Feller 
formula nor numerical calculations (22l.[23| are applicable. 

From the formulae (j4]) and ([5]) we can retrieve all ex- 
actly known cases enumerated in [T8l - [2ll | and give an 
unlimited number of new exact solutions g{a, (3;x); e. g . 
for a = 2/3, /3 = 1/3 (here ^ = 4, A: = 6, r = 1, see (Hi) 
and a; > 



2 1 
^l3'3'" 



c,(4,6,l) 



-2j73 



p , l,A(4,l+2i/3) 
5^6 I A(6,i-I-1) 



and for a; < 0, where ^ = 4, A; = 6, r = 3: 



(6) 



513,-3;- 



E 



c,(4,6,3) 

\x\l+23/3 



5^6 



1,A(4,1+2j73) 
A(6j + 1) 



66|x| 



(7) 



These two components neatly match at a; = 0. The coef- 
ficients Cj(4,6, 1), in Eq. (jG]) are equal to ^'"'^l^'^ 

J_ -5r(2/3) 7-3-"/" 
3tt ' 



37/2 



r(2/3) 



, 2r(2/3) n ^ n 14-3-"^ 
371- ' 37r ' r(2/3) 



9r(2/3) ' 

, and in Eq. ([7]) Cj(4, 6, 3) are equal 
respectively. This density 



is depicted in Fig. S) 

In Fig. [2] we present new distributions g(7 /5, x) for 
different (3, including symmetric g{7 /5,0; x), not explic- 
itly discussed in [20|. In Fig. 2] we display the PDFs 
g{a,l/3;x) for different a, including a ~ 2/3, see 
Eqs. ([6]) and ([7]), as well as g(l/3, 1/3; x) confined to 
—00 < X < 0. Our method confirms all the results 
for symmetric case /3 = obtained in related works 
[20| . in which unfortunately no graphical analysis 
was attempted. Our graphical representations warrant 
that it is a which grossly determines the global shape as 
well as their heavy tails. From Fig. [5] we also conclude 
that /3 has a much weaker influence on their shape as 
g{a = const, /3; x) for different /3 arc loosely evocative of 
each other. 

The precise knowledge of distribution functions 
g{a, P\ x) is a prerequisite to develop all the theories of 
anomalous diffusion based on the Fokker-Planck equa- 
tions l3Cl in conventional and fractional derivatives ver- 
sions [ItI . [30j The new solutions presented here offer a 



convenient starting point to carry out a systematic study 
of this approach, as they explicitly contain the parame- 
ter /? (usually set to /3 = —a in earlier attempts). Our 
solutions will also directly apply to analysis of hydrogen 
diffusion in the amorphous, high-temperature phase of 
PdssSiisHr.s [3] in which a = 1.54. In econophysical 
context [H the values a ~ 1.42 — 1.81 were observed, 
whereas in the values a = 1.64, 1.78 and 1.33 were at- 
tributed to the fits of statistics of 2000 Dow Jones Indus- 
trial Averages, to 1635 Boeing stock price returns and to 
the fluctuations of Yen/US$ exchange rate (1978-1991), 
respectively. 
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g(7/5, 3;x) 0.2 




FIG. 2. Comparison of asymmetric distributions g(a, P; x) for 
a = 7/5 and /3 = -3/5, -2/5, -1/5, 1/5, 2/5, respectively, 
starting from the right, 
case a = 7/5, /3 = 0. 



The dashed curve is the symmetric 



In conclusion, the presence of two natural parameters 
a and /3 in our solutions could permit a precise descrip- 
tion of experimental and statistical data characterized 
by Levy distributions. The index a governs the heavy 
tails whereas a and {3 adjust the position of distribution 
peaks. We expect that these solutions will be of use in 
the wide, boundary-crossing field of applications of Levy 
laws. 
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